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The paper describes perturbation techniques used for identification of rotating system 
dynamic characteristics. A comparison between two periodic frequency— swept perturbation 
methods applied in identification of fluid forces of rotating machines is presented. The description 
of the fluid force model identified by inputting circular periodic frequency— swept force is given. 
This model is based on the existence and strength of the circumferential flow, most often generated 
by the shaft rotation. The application of the fluid force model in rotor dynamic analysis is 
presented. It is shown that the rotor stability is an entire rotating system property. Some areas for 
further research are discussed. 


1. INTRODUCTION 


Perturbation methodology is well known and widely applied in Control Theory. The main 
objective of perturbation is identification of the svstem characteristics by exciting the system with 
a known input function and observing (measuring) its response. The comparison between the input 
and output of the system gives an insight into the system characteristics, i.e., its transfer function 
(Fig. 1). The "system" can be a mechanical structure, or any other system, such as biological, 
economic, etc. When little is known about the system characteristics, the system is often referred 
to as a "black box." When the mathematical structure of the system characteristics is known (or 
rather assumed approximately known, such as, for instance linearity of the n— th order), the system 
is referred to as a "grey box." The identification procedure provides the parameters to the assumed 
mathematical structure of the system model. The transition from a "black box" to "grey box" is 
performed by iterative adjustments to the system model. It results from a series of perturbation 
tests conducted for various operating conditions. 
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1.1 Perturbation Techniques in Rotating Machinery 

Rotating machines are particular subjects of perturbation testing. They represent a special 
class of "active" mechanical structures. The main motion of a rotating machine is rotation, driven 
by an external energy source. In order to fully understand the rotating machine process, and 
identify its operational characteristics, the perturbation tests should be conducted on a machine at 
its operational conditions, in particular when the rotor rotates at its operating speed (Fig. 2). The 
rotative energy provides, therefore, an additional input to the system. Most often this input is 
considered uncoupled from the perturbation input/output flow. The rotating energy, as an 
additional input to the system has, however, a very important influence on the system: its 

characteristics become rotative speed dependent. The perturbation input should be entirely 
independent from the rotative energy input. This technique is called "nonsynchronous" 
perturbation [1]. 

There are two particular cases of the perturbation technique applied to rotating machinery. 
In these cases only one input is used: (i) perturbation of the rotating machine at rest (like a 
"passive" structure), and (ii) synchronous perturbation. In the latter technique the rotative energy 
represents the only input to the system. In the synchronous perturbation method a controlled 
unbalance transfers rotational energy into an input force [2,3]. Both these particular perturbation 
techniques often provide important, but also very limited, information about the rotating machine 
characteristics. 


1.2 Characteristics of Rotating Machines 

Perturbation testing of rotating machines is used for identification of their mechanical 
characteristics associated with modes of vibration. Shaft lateral/ bending modes attract the highest 
attention. In particular, the lowest lateral/ bending modes are of the greatest concern because, due 
to the specific role of internal damping in rotating systems, they are usually characterized by the 
highest vibration amplitudes, thus creating conditions for rotor failures. 

In most rotating machines the shaft rotates in a fluid environment. Fluid involved in dy- 
namic motion becomes an important part of the system. Fluid/solid interaction causes the appear- 
ance of additional modes of vibration in the system. These particular modes have been identified 
by nonsynchronous perturbation testing [1]. Perturbation testing is used for identification of the 
rotating machine characteristics at various conditions such as under steady— state or variable load, 
unbalance, at different rotative speeds, including operational speed and above, at various 
temperatures, and other conditions. The information provided from the perturbation testing is used 
for model adjustments; it helps to predict stability of machine operation, provides tools for 
malfunction diagnostics, and assists in optimization of the machine performance. 

1.3 Perturbation Input 

Rotating machines can be perturbed by all classical input functions used in the identification 
of mechanical structures. Frequency swept periodic inputs [4—27], random inputs [28], and impulse 
inputs [29—35] became the most popular ones. Provided that the identified system is linear, all 
perturbation methods should lead to the same results. Each method has, however, its own strong 
and weak points which result in accuracy differences in specific applications. Two methods of 
nonsynchronous perturbation which have several advantages in comparison with other techniques 
are discussed below. 

1.4 Comparison of Two Frequency Swept Rotating Input Perturbation Techniques 
Used for Identification of Fluid Forces in Rotating Machines 

There are two main perturbation techniques of nonsynchronous one— mode testing used for 
identification of fluid force models in rotor/bearing, rotor/seal systems, as well as in fluid handling 
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Figure 3. - Nonsynchronous perturbation techniques: (a) perturbation by displacement (a 

stationary shaft and casing displacements are also used), (b) perturbation by force. 


machines (mainly pumps). The basic advantage of these methods: is the rotational .character of 
the input function in forward or reverse direction (the same or opposite to rotor rotation). This 
allows for independent identification of forward and backward modes of the rotating machine. 
These two perturbation techniques differ by the input /output functions (Fig. 3). For rotor circular 
orbits* they are as follows: 


DISPLACEMENT: Ae^* 


At(w) = 


HM 


FORCE: (F+jF,)e-* wt 


Refs. 7,11,12,18-26; 


FORCE: Fe^ . 

ftM = — 

DISPLACEMENT: Ae^ 4- ^ 





Refs. 4-6,8-10, 13-17; 


where A is rotor displacement amplitude, a is rotor response phase (also a = arctan (-Ft/F r )), ui is 
rotor perturbation (excitation) precessional frequency (usually varying from zero to some Wmax), t 

is time, j = -J—i, F r and Ft are radial and tangential forces respectively acting on the rotor. These 
forces are obtained by integrating the fluid pressure or measuring forces outside an elastically 
supported seal or bearing; F is the input force amplitude. 

In both cases the objective is identification of the "black box," i.e., the transfer function of 

the system. More precisely the sought functions are /c(w) = complex dynamic stiffness, or H(w) = 

— — = transfer function = complex dynamic compliance of the system. The most often obtained 
k(lo) 

result from either technique is the complex dynamic stiffness: 

F + jF 

k(u) = — when the input is displacement. 

k(u) = ^ e —ja when the input is force. (2) 



*A use of nondrcular forces/displacements or occurrences of noncircular response orbits/forces 
represent generalizations of the perturbation mehod, useful for higher eccentricities of the rotor [13]. 
For full identification of system parameters it requires a double number of measurements. 
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The complex dynamic stiffness components are as follows: 


Direct dynamic stiffness = DDS = = - C Q - S - 01 (3) 

F 

Quadrature dynamic stiffness = QDS = ■ (4) 

By limiting the input to a circular periodic function, both methodologies must yield exactly 
the same results, provided that the system is linear, and the instrumentation yields comparable 
signal— to— noise ratios. 

As is evident from Eqs. (3) and (4), the force— to— motion amplitude ratios F r /A and Ft/A 
correspond to the direct and quadrature dynamic stiffnesses, respectively, often, in some published 
papers, with opposite sign convention. The authors use the notion of dynamic stiffnesses, well 
established in Mechanics. For a system with one complex degree of freedom (one lateral mode of 
the symmetric rotor with no gyroscopic effect) the Direct dynamic stiffness is composed of the 
static (direct) stiffness (K), minus the effect of inertia and cross damping. By definition, the static 
stiffness K is positive; therefore, it appears in the direct dynamic stiffness with the plus sign. The 
quadrature dynamic stiffness contains the product of the radial (direct) damping (D) and frequency 
w, minus the constant (perturbation frequency independent) cross coupled stiffness. Since damping 
is positive, it is reasonable to present the corresponding straight line Do; versus perturbation 
frequency with a positive slope. The data from Refs. [21—23] illustrated in figures below have been 
presented in this unified format. The elements of the discussed above dynamic stiffnesses 
correspond to the standard linear bearing or seal models with isotropic properties (fluid inertia 
matrix with no cross terms, damping, and stiffness matrices skew symmetric). 

The results of perturbation testing of rotor/fluid systems are discussed in the next section. 


2. FLUID FORCE MODEL IDENTIFIED BY INPUTTING CIRCULAR 
PERIODIC FREQUENCY SWEPT FORCE 


As a result of over eight years of testing and identification research, the authors have 
proposed a fluid force model for lightly loaded bearings, seals, and fluid— handling machines, based 
on the significance of the steady circumferential flow [36]. The model was identified experimentally 
mainlv by using the frequency swept circular periodic perturbation force at the input of the 
rotor/bearing/seal systems. A similar fluid force model was previously developed theoretically, and 
has existed in a simplified version in the rotordynamic literature for at least 25 years (Bolotin [37], 
Black [38,39]). It has not, however, been fully exploited. Following the results of the perturbation 
testing, the most important improvement in the model introduced by the authors was the 
replacement of the constant "1/2," widely used in the fluid models, by the factor "A" (fluid 
circumferential average velocity ratio) as a function of fluid parameters and shaft eccentricity. This 
ratio has been identified for several cases of seals and bearings as a decreasing function of shaft 
eccentricity. 

The basic fluid circumferential average velocity ratio A is generated by the shaft rotation. 
Its value can, however, be strongly modified by forced fluid flow, such as preswirling and injections. 
Depending on the direction of the tangential preswirls and/or injections, A can be reduced when 
these directions oppose shaft rotation, or increased when the preswirls /injections are in the 
direction of shaft rotation. 
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For an originally concentric shaft rotating at a constant rotative speed U and perturbed by 
an unbalanced auxiliary system providing the frequency swept rotating force with frequency w in a 
range at least covering the interval OCuKfl for the "whirl" mode, the plot of direct dynamic 
stiffness of the fluid dynamic effects (shaft and other mechanical effects excluded), versus 
perturbation frequency uj approximately forms a parabola (Figs. 4 and 5), and thus can be described 
by three parameters. By using circular input forces and maintaining isotropy of the system (rotor 
gravity load balanced by centering springs), so that the response orbits are circular, the authors 
identified the fluid force direct dynamic stiffness as [1]: 

DDS = Ko-M^w-Afip (5) 

& 

where Ko is the fluid film static radial (direct) stiffness, is the fluid inertia effect, and A is the 

fluid circumferential average velocity ratio. Note that the "cross damping" (skew-symmetric terms 
in the damping matrix) appears in (5) as a function of the fluid inertia and Coriolis’ acceleration. 

The quadrature dynamic stiffness versus perturbation frequency usually forms a straight line 
(Figs. 6 and 7); thus, it can be described by two parameters. The authors have identified these 
parameters as 

QDS = D(w — Afi) (6) 

where D is the fluid film radial (direct) damping. The term Du> represents the passive effect of the 
shaft pushing the fluid, — DAQ is an active term transferring rotative energy into the fluid force 
pulling the shaft (tangential force). The quadrature dynamic stiffness (6) is the most representative 
and most important part of the fluid force model [40]. In spite of nonlinearities of both A and D as 
functions of shaft eccentricity, QDS appears as a distinctly straight line for all tested cases except 
for high shaft eccentricity [1,13]. This also means that there is practically no cross inertia effect in 
the fluid force. Note that in the identified fluid force model the direct and quadrature dynamic 
stiffnesses carry the common factor A. This means that the peak of DDS and zero of QDS occur at 
the same perturbation frequency (Figs. 5 and 7). The results obtained by other researchers show, 
however, some differences in frequency for these values (Figs. 4 and 6). This problem will be 
discussed in Section 7. 

Figures 4 to 7 illustrate the basic results obtained by various researchers, by using 
nonsynchronous perturbation techniques. 




Figure 4. - Direct Dynamic Stiffness versus 
perturbation frequency for the system with 

circular displacement input, and force out- Figure 5. - Direct Dynamic Stiffness versus per- 
put. 6 denotes a nondimensionalizing co- turbation frequency for the system with circu- 
efficient. Courtesy of H. Ohashi [22]. lar force input, and displacement output [10]. 
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The original response vectors versus perturbation frequency are presented in Figs. 8 and 9, in 
the form of Bode plots. Note that, in the technique which uses force input, the displacement 
response vector has a form characteristic of responses for a one— mode system to a periodic 
excitation with sweep frequency (Fig. 9). The occurrence of a resonance is obvious. In comparison 
with the response of a classical mechanical system, the response phase is, however, ahead of the 
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input force, and the phase sharply drops around 0° . This indicates the "quadrature" nature of this 
particular resonant phenomenon. The resonant frequency occurs at about AO, where 0 is the 
rotative speed. This means that AO represents one of the system natural frequencies , which is 
generated purely by the fluid/solid interaction. The system QDS becomes zero. This fluid— related 
quadrature resonance was documented by Stone & Underwood in 1947 [4], and again by Hull in 
1955 [5], but these excellent works were not immediately pursued. When the force is used as the 
system input, as opposed to a displacement input, the physical interpretation of the results becomes 
clear. The peak amplitude, which occurs when the quadrature dynamic stiffness equals zero, is of 
course limited by the total direct stiffness term in exactly the same fashion as the resonance at the 
zero of the direct stiffness, (yielding the classical stiffness— over— mass— related natural frequency) is 
limited by the total quadrature stiffness term (classical damping). The peak amplitudes may 
become very high, if the shaft rotative speed approaches the 1/A value of the rotor system 
"mechanical" natural frequency, corresponding to its first bending mode [16]. (Note that the 
responses include rotor mechanical effects.) In fact, instability occurs when both direct and 
quadrature dynamic stiffnesses zero at the same frequency. This subject is discussed in Section 5. 

By inputting a constant circular displacement and measuring the output in terms of forces, 
not only the accuracy of the results is lower, but also there is no clear physical interpretation of the 
results. The plot of force response amplitude versus perturbation frequency has an 
"anti— resonance" shape, a concave curve, a mirror image of displacement response amplitude versus 
frequency (Fig. 8). The phase is exactly the same in both techniques. It is obvious: by definition, 
•the phase represents the angle between the input vector and output vector, independent of the 
nature of the input and output functions. Note that the dip point of the response force amplitude 
occurs around one half of the rotative speed, i.e., the fluid circumferential average velocity ratio A is 
close to 1/2. The notion of the "force resonance," or rather "anti— resonance" is not well known in 
Mechanics. That is why by using the input force methodology, the authors have been more 
fortunate than those who used forces as outputs, to see direct physical interpretation of the results, 
namely that the value close to Afi is a system natural frequency. 

Note the advantage of the circular perturbation, as opposed to unilateral perturbation as in 
impulse testing. The results of forward perturbation (rotation and precession in the same direction) 
and backward perturbation (opposed direction) are significantly different (Fig. 9). The resonance 
occurs only for the forward perturbation, i.e., the "quadrature" natural frequency Afi of the system 
has the plus sign only (classical "direct" natural frequencies governed by stiffness and mass have + 
and - signs). 


3. FLUID FORCE MODEL IN DIFFERENTIAL FORM 


In more general terms, the identified fluid force model can be presented as follows. In the 
coordinate system which rotates at an angular velocity Afi the fluid force is assumed radial , with 
three components (Fig. 10): 

Fluid force in rotating coordinates = MfZ r + Dz r + K 0 z r (7) 

where z r (t)=x r (t)+jy r (t) are the shaft lateral displacements in rotating coordinates; Mf, D and K 0 
are the fluid inertia effect, damping and stiffness respectively. 

The coordinate transformation 

z f = z , where z(t) = x(t) + jy(t) , j=V I T (8) 
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from rotating z r (t) to stationary coordinates z(t) yields the fluid force: 


Fluid force in stationary coordinates = M^(z -2jAQz-A 2 Q 2 z) + D(z-jAQz) + K 0 z (9) 

The fluid force model, in differential form (9), can be directly used in rotor dynamic 
analyses. In particular, for shaft periodic circular motion with a constant amplitude: 

z = Ae^ (10) 

Fluid force for periodic motion = A-f-M^w-An) 2 + jD(u^-AQ) + K 0 }e^ = 


= Ae^ {DDS + jQDS} e Ae^ {s(u»)} 


( 11 ) 


where ac(o;) is the previously discussed complex dynamic stiffness. 


4. IDENTIFICATION OF THE NONLINEAR MODEL 


While determining the appropriate interpretation and use of the direct stiffness basic 
parabola, and the quadrature stiffness straight line, the authors are also investigating situations 
where the dynamic stiffness components diverge from parabolic and straight line shapes (Figs. 11 to 
13). These are cases of nonlinearities, and discontinuous shifts of various basic parameters. An 
example is shown below. 


•2 -I 0 12 3 



Figure 11. - Direct and Quadrature Dynamic 
Stiffnesses with nonlinear effects (perturbation 
by displacement). 6 denotes a nondimension- 
alizing coefficient. Courtesy of D. Childs [21]. 
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By inputting a force with sufficiently high amplitude F = mro> 2 (m,r are mass and radius of a 
controlled unbalance), the response becomes modified by system nonlinearities. The resulting 
dynamic stiffnesses diverge slightly from the parabolic/straight line shapes (Figs. 12 and 13). This 
occurs most noticeably in the range of perturbation frequency where the response amplitude is the 
largest (around a ; = Afl) (Fig. 9). By inputting forces with sequentially increasing amplitudes, F, 
for each consecutive test, the nonlinear functions in the fluid force model can be identified. This 
identification technique is graphically presented in Fig. 14. As a first step the DDS for a linear case 
(when F is small) is subtracted from the DDS for a nonlinear case. The result gives the stiffness 
nonlinear function ip versus perturbation frequency. Using this relationship, as well as the response 
amplitude versus perturbation frequency for high F , the plot of ip versus shaft eccentricity can be 



Figure 12. - Direct Dynamic Stiffnesses with Figure 13. - Quadrature Dynamic Stiffness with 
nonlinear effects (perturbation by force for nonlinear effects (perturbation by force) [17]. 

several cases of unbalance masses) [14]. m denotes mass of perturbation unbalance. 



Figure 14. - Identification of fluid film stiffness nonlinearity (perturbation by force). 
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obtained. Note that the response amplitude A represents a rotating radial eccentricity of the 
shaft. Assuming lateral symmetry, the rotating eccentricity corresponds to the unidirectional static 
eccentricity of the rotor. 

By a very similar method, the fluid film radial damping nonlinearity can be identified from 
the quadrature dynamic stiffness, provided that the nonlinearity of A is relatively small [17]. 

Note that the nonlinear effects in the fluid film may be caused and/or influenced by fluid 
injections and preswirls [21], as well as cavitations [26]. 


5. ROTOR STABILITY: CRITERION AND MARGINS OF STABILITY 


The rotor/fluid system dynamic stiffness components as functions of perturbation frequency 
can be used to predict the rotor stability and the stability margin. 

In very general terms, the threshold of stability occurs when the "zeros" of the direct and 
quadrature dynamic stiffnesses occur at the same perturbation/precession frequency. In terms of 
frequency, the stability margin is the lowest difference between the zeros of the direct and 
quadrature dynamic stiffnesses (Fig. 15). 

In actual applications, both components of the fluid force dynamic stiffness become modified 
by the rotor system parameters. An example is discussed below. 

Consider an isotropic rotor supported in one rigid and one fluid— lubricated bearing (Fig. 16). 
The rotor rotates at a constant speed Q and is perturbed by a nonsynchronous periodic circular 

forward (in the direction of rotation) sweep frequency force Fe^. The model of the system is as 
follows: 



Figure 15. - Rotor stability prediction using Figure 16. - Rotor/bearing system model, 

dynamic stiffness graphs. 
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Mz + Mj(z -2jAt!z -AWz) + D s z + D(z -jAfiz) + Kz + K 0 z = fJ" 1 


( 12 ) 


where M, D s , K are rotor generalized (modal) mass, damping, and stiffness, for the first lateral 
mode respectively. The rotor system response is circular: 


z = (13) 

where the amplitude A and phase a can be calculated when Eq. (13) is substituted in Eq. (12): 

Ae ja [-M£ t ;2_M f (a^Afi)2+D s ja;+Dj(a>-An)+K+Ko] = F (14) 


From Eq. (14) the amplitude A and phase a of the forced response (13) are obtained as 
follows: 


/ [K+K 0 -Mw 5 HVI f (a>-Ai2)2]2 + [D g OM-D(^Afl)] 2 


(15) 


D UH-Dfw- XCl) 

a = arctan Mw 2 +Mf ( ) 2 -k- K 0 

The whirl resonance (Fig. 9) occurs when the quadrature part 
is equal to zero, i.e., when 

Afi _ , 
w l+D g /D = W res 

The resonant amplitude and phase are 


(16) 

of the denominator of Eq. (15) 

(17) 


A res |K+K 0 — (M+M f D2/D 2 )A 2 n 2 /(l+D g /D)2| ~ | K+Ko-^lA 2 i2^| » “res 

The actual peak amplitude depends, however, on whether the force amplitude F is constant 
or dependent on frequency u. 

Eq. (16) indicates the leading-phase phenomenon for low frequency. At a; = 0 the phase is 
a = arctan > 0 

IV-fA-O 

For AQ > 0 the response is, therefore, ahead of the exciting force. 

Eqs. (15) and (16) also show clearly that the whirl resonance occurs only if Q and u have the 
same rotational directions (perturbation forward). If u < 0 (perturbation backward), then the 
quadrature resonance does not take place (compare with experimental results, Fig. 9). 

The system complex dynamic stiffness is obtained from Eq. (14) as the ratio F/(Ae^“). The 
dynamic stiffness components are 
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DDS = -Mo;2-M f (a^Afi)2 + K + K 0 

QDS = D w + Dfcth-An) 
s 

The system dynamic stiffness is, therefore, a combination of the rotor and fluid— related terms. 
The eigenvalue problem for Eq. (12) provides the following eigenvalues: 

Si, 2, 3, 4 = 2(M+M ) (-D-© 8 +2jAfiM f ) ± -j: [J -P+VPTR? ± ij j 

f 


(18) 

(19) 


where 
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K+K o— M^O 2 


M+M_ 


D+D * 
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(M+M f )2 M+M f 


Afi 


The stability criterion, i.e., the requirement that all eigenvalues have non— positive real parts 
is as follows: 


D+D 


-2(M+4 + ^'^^ < 0 

which yields the rotative speed threshold of stability: 
D 


n < 


1+ 


u 

a 

D 


K+K p 

M+M.D2 /D2 
I s ' 


( 20 ) 


The natural frequency at the threshold of stability is 
part of the eigenvalue). 

Inequality (20) can also be solved in another form: 


m+m K dW ( this is the ima e inar y 

f s' 


An 


l+M/M 


K+Ko 

M+M. 


A2Q2 


MM, 


An 


An 


K+K 0 MM f 


i- y f .».T^ U A 20 2 

(M + M f )2 - 1+D s /D - 1+M/M f + / M+M f (M+M f )2 


The direct dynamic stiffness (18) is equal to zero when 


( 21 ) 


uj = w, 


Dl,2 l+M/M 


/S^-A2n2 


MM, 


M+M, 


(M+M f ) 2 


( 22 ) 


The values u are called roots (zeros) of the direct dynamic stiffness. The quadrature dynamic 

U lj it 

stiffness (19) is equal to zero when 


u= UJ n = 


An 


Q " 1+D S /D 


(23) 
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where u is the root (zero) of the quadrature dynamic stiffness. The stability criterion (21) can, 
therefore, be written as 


LO < < UL 

D 1 - Q “ D2 


which means, for rotor stability, the root of the system quadrature dynamic stiffness should occur 
between the two roots of the direct dynamic stiffness. 


roots: 


The margin of stability , in terms of frequency, is defined as a minimum distance between 


(24) 

This margin of stability is "horizontal," thus expressed in terms of frequency (Fig. 15). An 
additional "vertical" margin of stability can also be defined. It expresses the direct and quadrature 
dynamic stiffness closeness to the threshold of stability (Fig. 15): 

min [DDS(at w Q ), QDS(at o> Dl ), QDS(at ^ D2 )] (25) 


The stability margin (25) indicates which parameter of the system is most likely responsible 
for the instability. 

The threshold of stability occurs when both direct and quadrature dynamic stiffness have the 
same root, i.e., either w = u or u — u> . It means that at the precession frequency when these 

l) 1 VJ y Is it 

equalities occur, the total complex dynamic stiffness of the system equals zero (see Eq. (15) in 
which the denominator is zero at this case). 

When the threshold of stability is exceeded, the vibration amplitude increases, most often 
ending up in a limit cycle self-excited vibrations, known as "whirl" and "whip," which are governed 
by the system nonlinearities [40, 41]. 

The fluid force model presented in the dynamic stiffness form provides clear suggestions for 
stability improvements of rotor/bearing/seal systems. Higher stability margins result when 

• A decreases (provided that DM>D s Mf, which is almost always satisfied) 

• Ko increases 

• D decreases (most often is, however, insensitive to stability) 

• Mf decreases 

• K and D s increase 

• M decreases 

A higher stability margin results also when the rotative speed fi decreases. Since variations 
of fi can involve other rotor dynamic effects, such as balance resonance which may interact with 
fluid effects [40], 0 is not included in the above list. 

The most effective parameters for rotor stability control are a decrease of the fluid 
circumferential velocity ratio A and an increase of the fluid film radial stiffness Ko. The first is 
widely known in "anti— swirl" applications [41—51], The second, i.e., the increase of Ko results 
directly from an increase of shaft eccentricity and/or increase of the fluid pressure. The latter 
conclusion applies to externally pressurized ("hydrostatic") bearing and seal designs, which are 
widely known for their stability features. The most surprising result concerns fluid radial damping 
D . In most cases, an increase of D has no influence on stability, or has even a minor negative 
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effect. (This effect can be compared to destabilizing action of rotor structural/ internal friction at 
high rotative speeds.) 

An interesting result is obtained if the classical Half Power Bandwidth method is applied to 
the whirl resonance. The classical amplification factor Q may be calculated as 

q _ k^res 
^ U)2 - U3\ 

where a; res is given by Eq. (17) and W4, w 2 are frequencies corresponding to the phases equal to 

±45° or amplitudes at A (3 dB level). For the whirl resonance the amplification becomes 

res' 

(Fig. 9): 

Q M S 1 ( 26 ) 

D V(K+ K o)/M /Afi - 1 

The first fraction of Eq. (26) is the classical relationship of half the inverse of the damping 
factor. The second fraction indicates the closeness of the direct resonance (zero of DDS) to 
quadrature resonance. (Shaft damping and fluid inertia effect were omitted in Eq. (26) for clarity.) 
The whirl resonant amplitude becomes significantly magnified if the rotative speed is close to 1/A 
times the rotor lateral mode natural frequency. This effect has often been observed experimentally 
[16], 


Finally, the dynamic stiffness format can be used to successfully predict rotor instability in 
case of nonlinear distortions of the basic parabola/straight lines. The general idea always holds 
true: Instability occurs when both DDS and QDS have zeros at the same perturbation/precession 
frequency. Nonlinear effects may cause dramatic changes in the shape of DDS and/or QDS, 
including multiple zero points. Theoretical studies [52] proved, for instance, that the low frequency 
rotating stall instability can be predicted using the method discussed above. 

Note that the relationships discussed in this chapter hold true for the considered example of 
one mode isotropic rotor/bearing or rotor/seal system. When more modes are taken into account, 
the relationships become more complex. 


6. APPLICATIONS OF THE FLUID FORCE MODEL 


By implementing the model identified by the authors into rotating machine dynamic 
equations, several significant and valuable results were obtained. The results of analysis conducted 
for one, two, and three mode models of rotors are very encouraging [36,40,41,42,53,54]. The model 
allowed for explanation of new phenomena, such as the recently observed second mode whirl [53] 
and second mode whip [54]. The generalization from two and three modes to "n" modes is obvious 
[53]. All swirl preventive measures (such as anti— swirl, preswirl, swirl brakes, fluid injections) can 
also be easily included in the models [41,42], and thus would give designers very good insight into 
the expected efficiency of such devices for stability of rotating machines. 


7. FLUID FORCE MODEL ADJUSTMENTS AND GENERALIZATIONS 


The model of fluid force (9) identified by the authors, using force input perturbation testing, 
is certainly open to further developments. One adjustment, namely the addition of nonlinear terms, 
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was discussed in Section 4. There exist also some other possible modifications to the model. These 
include: 

(i) Tangential components. In the original form (7), the fluid force was assumed to have 
radial components only. There is a high probability, especially when the shaft eccentricity 
increases, that the force components are not exactly radial, but also have a tangential part. For 
instance the stiffness force in (7) may require an addition: the term Kqz t should be replaced by 
(K 0 -jK t )z r where Kt is tangential stiffness. 

(ii) Two fluid circumferential velocity ratios. In the model (7), it was assumed that all 
components of the fluid force rotate at the same angular velocity AQ. It might be necessary to 
introduce two different ratios Ai, A 2 , one for the fluid inertia force, the other for the damping force. 
Thus, the adjusted model (9) will contain the terms: 

M^z- 2jA jOz— Ai$l 2 z) + D(z— jA 3 Dz) . 

where Aifi and A 2 U are angular velocities at which inertia and damping forces rotate 
correspondingly. Two different A ratios seem to be required in order to adequately identify the 
fluid force model in pumps (see Figs. 4 and 6 [22]). 

< (iii) Higher order terms. Some results of perturbation testing [26] suggest that the fluid 
force contains terms of order higher than two. In particular a third order term (jerk) should be 
included. With an assumed additional value of the corresponding fluid circumferential velocity 
ratio, A3, the model (9). might, therefore, be completed by the following term: 

J(z-3jA 3 f2z -3A 3 fi 2 z + jAjjfi 3 z) 

where J can be considered as a complex parameter, containing radial and tangential components. 

(iv) Nonsymmetric fluid force. Most results of the perturbation testing have been obtained 
for symmetric cases: perturbation around the shaft centered position. When the shaft is statically 
displaced to higher eccentricity, the flow conditions are certainly modified. The fluid force then 
contains terms of nonsymmetric nature [13,14]. In particular, the nonlinear terms identified for 
symmetric cases as simple functions of radial eccentricity, may require adjustments. More research 
is required in this area. 


8. CONCLUDING REMARKS 


This paper discusses the use of nonsynchronous swept frequency harmonic rotating input 
perturbation techniques for identification of the fluid force model in lightly loaded 
rotor/bearing/seal systems. Two basic methods differ by a choice of input functions: either as a 
orce, as used by the authors and a few other researchers [1,2,4—6,8—10,13—17], or as a displacement 
7,11,12,18—26]. The measured output responses are displacements or forces respectively. While 
)oth methods provide the system dynamic stiffness characteristics, the main difference between 
these two methods consists in data generation, acquisition, and processing; thus they provide 
different levels of signal— to— noise ratios. A significant contributor to noise is input /output phase 
readings. The first method provides a higher accuracy (around ±0.3%) than the second method, 
especially in phase measurements. An estimated error of the first method final results is lower than 
±3%. An application of both methods to the same system (if possible) might ultimately provide the 
best identification data. 
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The results, in terms of the system dynamic stiffness components obtained from both 
methods, are very similar, as should be expected. Based on results of hundreds of tests obtained 
over the past eight years, the authors have provided a consistent mechanical interpretation of the 
results, and proposed an improved model of fluid forces in lightly loaded bearings and seals [36]. 
The model includes the fluid circumferential average velocity ratio as a function of shaft 
eccentricity, instead of a constant value 1/2 as in the classical models. (In some recent analyses the 
assumption about the constant 1/2 ratio starts gradually fading away.) The bearing/seal 
coefficients can easily be obtained from the proposed model as particular cases. The application of 
the improved fluid force model in rotor/bearing/seal systems yields results which stand in a very 
good agreement with observed and documented dynamic phenomena, such as whirl, whip, and 
higher mode whirl/whip vibrations [36,40,41,53,54]. The model allows for an improved prediction 
of the rotor stability thresholds and stability margins, as well as quantifying stability control 
measures. Some possible further adjustments of the fluid force model are also discussed. The latter 
would require more, well coordinated, experimental research. 

The database resulting from different types of perturbation methodologies applied to 
rotor/bearing/seal/pump systems is now quite rich, including data from CalTech, Mitsubishi, 
Sulzer Brothers Ltd., Universities of Kobe, Tokyo, Osaka, Kaiserslautern, Texas A&M, and Case 
Western. It will be very useful to reduce the acquired experimental, as well as 
analytical/numerical, data in the form proposed in [55]. It will then become clearer what 
adjustments the model would require. It will also indicate directions for further research. 
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NOTATION 


A 

D,K 0 ,M f 
D S ,K,M 
DDS, QDS 
F 

F r , F t 

s 

t 

z = x+jy 
z r = x r +jy r 
a 


Amplitude of rotor vibration 

Fluid radial damping, stiffness, and inertia effect respectively 
Rotor generalized (modal) damping, stiffness, and mass respectively 
Direct and quadrature dynamic stiffness components respectively 
Force amplitude 

Radial and tangential force components 

Eigenvalue 

Time 

Rotor displacement in stationary coordinates 
Rotor displacement in rotating coordinates 
Force/displacement phase 


K 

A 

U) 

n 


Complex dynamic stiffness 
Fluid circumferential average velocity ratio 
Perturbation/precession frequency 
Rotative speed 
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